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ABSTRACT 

The 0{a^^) and 0{aas^^) corrections to the partial decay width T{Z — > 6 + 6) are 
computed using the current algebra formulation of radiative corrections. This framework 
allows one to easily enforce the relevant Ward identity that greatly simplifies the calculations. 
As a result, the one- loop 0(a— ^) contribution is computed through the investigation of only 
two convergent diagrams. The computation of the QCD corrections to the one-loop 0(a-^) 
term involves fewer diagrams than in the standard approach. In particular, the number of 
infrared divergent contribution is reduced. The calculation is performed in the dimensional 
regularization scheme and no term more divergent than is found. Our result confirms 
the screening of the one- loop top mass effect recently found by Fleischer et al.. 
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1 Introduction 



It is well known that the electroweak radiative parameters [Q] and Ar are affected by 
virtual top exchanges in the self-energies of the vector bosons through terms that depend 
quadratically on mt, the mass of the top quark. Moreover, due to the isodoublet nature of 
the top, corrections are also found in the Z66-vertex. 

The status of the theoretical calculations of these corrections is quite advanced. At 
the moment we have available for the vacuum polarization functions, beside the one-loop 
calculation the complete perturbative 0{aas) QCD contribution studies of the ti 

2 _ 

threshold effects F] and the leading 0((a^-)^) term [0,|8|. Concerning the Zbb-vertex, the 
one-loop calculation was performed by several groups a few years ago [p|-pT[|. Recently, the 

2 2 

leading 0{aas-^) ||l2i and the 0((a^^)^) terms H have been computed. 

The calculation of the two-loop corrections has proven to be a difficult task. Indeed, 
apart from the perturbative QCD contribution to the vacuum polarization, the two-loop 
contributions have been evaluated in some approximation, either neglecting all the masses and 
momenta but the top mass [[^,p^ or retaining besides the top only the higgs mass [^]. However, 
these approximations are sufficient to derive the most interesting part of the corrections, 
namely the leading mf contribution. 

It is the aim of this paper to apply techniques used in the current algebra formulation of 

2 2 

radiative corrections to derive the leading 0{a^-) and 0{aas^-) corrections in the 
partial decay width T{Z b + b). As will be seen in the following, these techniques greatly 
simplify the calculations. In the most interesting case, namely the two- loop 0(aas-^) 
correction, in addition to facing fewer contributions than in the standard calculation |]l2| one 
deals with less divergent diagrams . Furthermore, the treatment of the infrared divergent (IR) 
contribution is greatly simplified. Our result shows that the QCD corrections have opposite 
sign to that of the one-loop contribution, leading then an increase in the mt upper bound. 

^For an updated analysis of the QCD effects in the electrowealc corrections see Ref. and references 
therein. 
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We find agreement with the recent calculation of the asG^m^ term by Fleischer et ai. 

The paper is organized as follows. In Section 2 the current algebra derivation of the 
0(0—7") one-loop term is presented. Section 3 is devoted to the calculation of the QCD 
corrections to this term. In Section 4 we discuss the results. The Appendices give some 
details about the QCD calculation. In Appendix A we show that the formulae used for 
the four and five-point correlation functions (Cfr. Eqs.(18ajl^)) contain correctly the two- 



loop wave function renormalization of the external legs. Appendix B lists the individual 
contributions of the various diagrams. 

2 

2 One-loop 0{a^^) correction 



In this section we derive the leading one-loop 0(a^-) term in the vertex using two and 
three-point correlation functions. This allows us to set the framework for our subsequent 
discussion of the QCD corrections and to derive the basic Ward identity that enters into the 
calculation. 

In order to fix our notation we write the part of the Standard Model (SM) Lagrangian 



density that describes the interaction of the W, Z and unphysical scalars with fermions as |13] 



Cm = -^{WUI^ + h.c.) - ^-Z^Ji'z - TT^ [^2^2 + V2i<^^S + h.c.) 
V2 c zniw L 



where g is the SU{2) coupling, niw stands for the mass of the W boson, c is an abbreviation 
for cos 9w, Jz ^^'^ fermionic currents coupled to Z and W respectively, is the 

field that creates a meson, $2 and $ the unphysical counterparts associated with the Z 
and W and 

S2 =2 df,J^ = -ii) mO CsTsV' (2a) 
S = -id^J>^ = ijT^I,. (2b) 
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In Eqs.(pa|j2bD ip represents the column vector ip = [t, b)'^ , mP , C3 and F are the 2x2 matrices 



m 



ml 




(2c) 





-1 



(2d) 



r 



(2e) 



a+ = 



1+75 



and the superscript on nit refers to the bare mass. As it is evident from 



Eqs.( p^ , p^ ) we are considering only the third generation and taking the bottom quark as 
massless. 

2 

The one- loop vertex diagrams contributing to the 0(a^-) correction are depicted sym- 
bolically in Fig.l. The circles represent the sum of diagrams in which the ends of the <I>^ 
propagators are attached in all possible ways to the external bb quarks. 

The amplitude corresponding to Fig. 1(a) can be expressed as E 



= - i - lim T^{q, p, p') 
c q^q 



(3a) 



with 



-I r 



(Pk 



2 2 ml, J (27r)'^ fc2 _ ^ 



(Tye 



-iq-y 



X / d'^xe'^-^ < p'p\T* \j^{y){S^x)S{0) + h.c.) 



0>, 



(3b) 



where n is the dimension of space-time, fi is the 't Hooft mass scale, T* is the covariant 
time-ordered product, p and p' the momenta of the b and b quark respectively and q = p + p' 
the momentum carried by the Z. As the 0{a^-) term should be gauge invariant by itself 
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we find it convenient to carry out the calculation in the 't Hooft-Feynman gauge in which 
the propagator of the ^ field has the form i{k'^ — ml 



As shown in Ref. |13, 14|, the limiting procedure in Eq.(^) affects the insertion of the $ fields 
on the external lines in such a way that contains not only the proper vertex correction 
but also the contributions from the wave function renormalization of the b quarks. 
To trigger a Ward identity we contract with qa obtaining 

where 2?$ is an expression analogous to with the replacement ~^ ~'''d^J^ and 

V^{k) = --^^^^ frxe"'-^ <p'p\T*{S\x)S{0)-SH0)S{x))\0>, (5) 
2 J 

with = 1 — c^. In deriving Eq.(||) we have used the commutation relation 

6ixo - yo) \4{x), S{y)] = -^^S{x)6^{x - y). (6) 



Differentiating Eq.(0) with respect to q^j, one obtains 



It is easy to show following, for example, the discussion in Section VII of Ref. ||T^ that the 

2 

first term in the r.h.s. of Eq.(0) is 0(a) instead of 0(a^^). Terms of this kind will be 
considered throughout all our calculations as subleading and therefore neglected. After a 
partial integration in the last term in Eq.(|^) we can write for 



X I (Txe'^^-^^'' <p'p\T*{S\x)S{Q)-S\Q)S{x))\Q> . (8) 
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As we are interested in the leading rrit dependent contribution, we can set q = q = in the 
second term of Eq. (P) and obtain 



d% * 2 ml, 



where the subscript L.T. reminds us that we are considering only the leading term. 
Turning our attention to Fig. 1(b) we write the amplitude as 



(9) 



lA^ — — i—Ug, 

c 



9 



2m% 



(27r)" (A;2 - ml) [{k + qf - ml] 
J (Pxe'^-'' < p'p\T*{S\x)S{0))\0> . 



(10a) 



(10b) 



It is immediate to see that the leading term in Uj^, obtained by putting g = in Eq. (|10b|) , 

2 

cancels exactly the second term in the r.h.s. of Eq.(^). We conclude that the leading 0{a-^^) 



correction in the Zbb-vevtex is given by 



d 



'8c ml dq^ 

J,k-x 



;9 d 



(F-k 



q=q 



or, recalling Eq.HZ 



(27r)" ;u"-4 A;2 - ml 



rye- 



iq-y 



> 



(11) 



q=q 



To express the r.h.s. of Eq.(|ll|) in terms of Feynman diagrams we use Wick's theorem. 
We now observe that a non-zero contribution is obtained only when 5*2 is acting on a top line. 
Fig. 2 represents schematically the fermion lines in Eq.(^) for q ^ q before the k integration 
is performed. The dotted line stands for the q momentum absorbed by the 5*2 operator while 
the dashed lines indicate the momenta emitted or absorbed by the S and operators. For 
example, in Fig.2(a) q is absorbed by S2{y) while k is emitted by S'^{x). Similarly in Fig. 2(b) 
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k is emitted by S{x) and q is absorbed by 5*2 (y). 

Using an anticommuting 75 the contributions of Fig. 2(a) and Fig. 2(b) to Eq.([Tl]) are 

mf d \ f d'^k u{p')^a-v{p) 

2(a — 



2(b) 



8c m'^ dq^ 17 (2vr)4 [fc^ — m^] [{k — p + qY — mf] [{k — pY — m^] J 

mf d \ f d'^k u{p')^a-v{p) 
8^^df^ \J (27r)4 [A;2 - m^] [{k + p' - qf - mj] [{k + p'^ - mf\ 



(12a) 

q=q 

(12b) 

q=q 



where the l.h.s. in the above equations indicates the appropriate diagram in Fig. 2 and 
a_ = ^^2^ ■ It is easy to see that the differentiation with respect to of the denomina- 
tors in Eqs.(|l2|) gives a zero contribution after the limit g — > g is taken. Differentiating the 
numerators and then putting p' = p = m-^ = one obtains 

{V'^+U^^)l.t. = -i-^-f^^uip'h'aMp)- (13) 
2cl67r"^2m^ 

Eq.([l^) gives the contribution of the leading 0(a^-) vertex correction. In the on-shell 
scheme where = sin^ 6w is defined according to 

2 

2 

sin 6w = 1 



and the partial width T{Z f + /), where / is a massless fermion, can be written as [15| 
r(Z ^ / + /-) = {1 - ^ifQfs'Re kf + SQ)s^\kf\'] (14) 



we see that the term ( p!3| ) introduces a modification in the electroweak form factors pf,fe 
and /cf, equal to — 4xj and +2xt respectively, where xt = {G m^) / {\/2 Stt'^) . In Eq.([l^ Nc 
refers to the color factor, /| and Qj stand for the /s and charge quantum number of the 



fermion and Spirr = 3xj. The counterpart of pbb and k^ in the MS scheme, namely the 
quantity pf,b and kf, |15|, will be accordingly modified by the factors — Q!/(47rs^)(m^/m^.) and 
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a/(47rs^)m^/(2 77?.^,,), where a = e^(m2)/(47r) and = sin^ Owimz) are the MS couphng 
and weak angle evaluated at fj, = ruz- 



2 

3 Two- loop O(aas^) correction 

2 2 

We proceed to compute the QCD corrections to the leading 0{a^^) term, i.e. the Oiaas^^) 
contribution. We begin by noticing that the wave function renormahzation of the external Z 

2 

has no contribution 0(aa<j^^), therefore all the QCD corrections are obtained by adding a 
gluon, in all possible ways, to the diagrams of Fig.l. The result is shown in Fig. 3. The circles 
represent now the sum of diagrams of various kind, i) Diagrams where both the ends of the 
$ as well as the gluon propagator are attached to the external 6, b lines. These are reducible 
diagrams contributing to the wave function renormahzation of the external fermions. ii) 
Diagrams where either the $ or the gluon propagator has both ends attached to an external 
line while the other propagator is acting internally, iii) Diagrams where each propagator has 
one end attached to an external fermion while the other end is attached internally. 
Fig.3(a) can be expressed as a five-point correlation function or 



X y rx2e''''-''' J rxse-"''--^ <p'p\T* [j^{y)J^{x2) Js xixs) {8^x1) S{0) + h.c.) 



|0>, 
(15b) 



where gs and Jg are the SU{3)c coupling and current respectively and we have suppressed the 
color indices. Analogously to the one-loop vertex function V^, V^^ contains not only the one- 
particle irreducible (IPI) two- loop corrections to the proper vertex but also the contribution 
of the wave function renormahzation of the external lines. Concerning the latter, as shown in 
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Appendix A, the limiting procedure in Eq.(15a) does correctly take into account the two-loop 
IPI contribution of it but there is a mismatch in the numerical coefficient of the reducible 
part. To cure it we replace Eq.(15a) by 



V^^ = -z^ lim(T|:^-n«r^) (16) 

where 

m) = ^\ r(2 - -) (2 - n) C da (1 - a) \\^{\ - a)l-2+"/2 (17) 
3 (167r2)"/4 ^ 2^ ^ Jo 

is the QCD one- loop wave function renormalization constant in n dimensions and A is a 
fictitious gluon mass introduced to regularize the IR divergencies. 

The four-point correlation function depicted in Fig. 3(b), l^!^g, presents an analogous over- 
counting in the field renormalization of the external legs. Analogously to Eq.(|T^) the correct 
factor can be obtained by subtracting a term proportional to the one-loop two-point correla- 
tion function U^. We write the corrected amplitude as 

Z?^^ = -i^(<^-n«Z^^) (18a) 

7> = 2. 2 _ 2. f d^h K f rf"fc2 1 

*3 2ml^'^'' ' ' J (27r)" {kj - ml) [{h + q)^ - ml] J (27r)" 



|0 > . 

(18b) 



The terms 11^^' in Eq.(|T6|) and Wg^Uj^ in Eq.( |18a| ) can be understood as follows. In the 
one-loop vertex functions and U^, we can consider the bra < p'p \ as dressed with respect 
to the strong interactions. Expansion up to 0{as) gives the above two contributions. 

In order to reduce Eq.( |15b| ) to a more tractable expression we apply the same procedure 
developed in Section 2, namely we contract T^g — Wg^T^ with qa and then we differentiate 
with respect to q^. The important point to notice is that commutes with J^. Therefore, 
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going through the same steps as in Section 2, we reach the two-loop counterpart of Eq.(0) or 



in. 



d 



2/ 2 



(1^X2 e 



w 
ik2-X2 



(27r)" {kj - m2,)2 J (27r)" /i"-^ j^: 



rx^e-"''--^ <p'p\T* \j^{x2)Jsxix3)S\xi)S{0) 



0>, 



(19) 



where V^g is again an expression obtained by T^^ with the replacement — > —id^J^. The 



second term in the r.h.s. of Eq.(19) cancels the leading contribution of Eq.( 18b| ) and we are 
left with 



(20) 



q=q 



The diagrams contributing to are obtained by adding a virtual gluon, in all possible 
ways, to the graphs of Fig. 2. The types of diagrams are shown in Fig. 4. To every graph in 
the figure there corresponds two diagrams. In the first one the momentum ki is emitted by 
S^{x) while in the second it is emitted by S{x). 

To get an ultraviolet (UV) finite answer we have to add to Eq.(^) the counterterm 
contributions, that we indicate as c.t.. The only counterterm at our disposal comes from 
the bare top mass. Taking as renormalized quantity the zero of the real part of the inverse 
propagator, the so-called "on-shell" (OS) mass, we have for the counterterm 6mt 



Suit 



, ,rmtY{2--) I da[2(l-a)+nl(m2a2)-2+"/2. 

(167r2)"/4 * ^ 2^ Jo ^ ^ ' \y t ) 



(21) 



In Fig. 5 the types of diagrams contributing to c.t. are depicted. The cross represents the 
insertion of 5mt. Figs.5(a,b,c) are due to the presence of in the operators S2 and S 
(Cfr. Eqs.( p^j2e[ )). Similarly to Fig. 4, every graph in Fig. 5 represents two contributions with 
different momentum insertions. 
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Evaluation of the contributions of Fig.4 and Fig.5 (Cfr. Appendix B) gives 



where Cj? is the eigenvalue of the quadratic casimir operator for the fundamental represen- 
tation of SU{N), namely (iV^ _ l)/{2N), [Cp = 4/3 for SC/(3)c), 

ci = - ^ + (23a) 
C2 = - 6 (23b) 



and 



/"^ , ln(l - x) vr^ , , 

^(2) = - / dx^ = —. (23c) 



In Eq.( |23a| ) B^. represents the contribution of the IR divergent diagrams 4(a), 4(b), 4(c) and 
once the UV pole has been subtracted. Explicitly 



q=q 



Br = -\n^'-^-^\n'-^-Z+'—. (24) 



A2 A2 ^2 

Y-31n^-3 + — 

2 

Comparing Eq.(p^) with Eq.([l3|) we see that the QCD corrections modify the leading 0(a-^) 
term in the Z66-vertex by 

^ (ci + C2^(2)) Cp. 

2 _ 

The inclusion of the 0(aas— j-) corrections in the decay width T[Z b + b) can be 
performed in the following way. We write the vertex, V^, as 

= -^^^M^±^ (25) 
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where 



O =Oo + aie + ais + 02 



b =bo + bie + bis + &2 



(26a) 
(26b) 



with 



ao 
aie 
ais 

a2 



■- h -s 

2 3 

16tt^ 4 m?. 



2 . 



as 
in 
as 
in 



aie(ci + C2e(2))CF ; 



The decay width can be written as 



On =- 
2 

ble = — aie 

bis =^boCiCF 
47r 

b2 = - 02- 



T{Z ^b + b) = ^m^{a^ + b^). 
4o7r 



(26c) 
(26d) 
(26e) 
(26f) 



(27) 



Expanding the squares in Eq.(^), keeping terms up to the order we are considering, one 
obtains 



T{Z ^ b+b) 



487r' 



[al + bl){l + ^ciCf) + 2(ao - 6o)aie U + ^(2ci + C2m)CF 



(28) 



It is clearly understood that in Eq.(2S) the factor present in ci is evaluated at q 



A few remarks about the above result are now in order, i) All the calculation has been 
performed analytically. The integration over the Feynman parameters has been checked, 
always analytically, with the algebraic manipulation program MAPLE ii) In every 

diagram in Fig. 4 and Fig. 5 but 4(a) we have neglected all the momenta and masses compared 
to the top mass. In 4(a) we were forced to keep in addition the momentum transfer q due 



11 



to the presence of IR divergent terms, iii) We found it advantageous to regularize the IR 
divergencies by giving a small mass to the gluon instead of using dimensional regularization. 
In this way we did not introduce any term (e = n — 4) in the calculation keeping the 
n-dimensional algebra simpler. 



Eq.(|28|) is IR divergent in the ci coefficient (Cfr. Eqs.(23a) and (|24|)). As is well known, to 
eliminate the IR divergencies it is necessary to take into account the QCD bremsstrahlung. 
In particular, to be consistent with the order of our calculation, we need to consider the 
bremsstrahlung with the inclusion of the leading 0(a-^^) term. The discussion of this cor- 



rection can be performed on the same footing as the 0(aas-^) term. In fact, the cancellation 
between Eq.( [L8bD and the second term in Eq.(]l9|) will work even in the presence of a single 

2 

Jg operator. The result is that the 0(a^^) term in the QCD bremsstrahlung is obtained 
by adding a real gluon in all possible ways to the diagrams in Fig. 2. We note that a gluon 
emitted by a virtual fermion line gives rise to a correction that does not contribute to the 
leading term. It is then easy to show that the 0(a-^^) correction can be absorbed in a 
redefinition of the vector and axial coupling of the Z. Therefore we can write for the rate 



3 

-^^ + 2 



Cf 



(29) 



with 



d =aQ + oie 
b =bo + bie. 



(30a) 
(30b) 



Expanding the squares in Eq. (^9)) , retaining terms up to the order we are interested, we have 



T{Z ^b + b + g) = ^m^ 



(al + 6o + 2(ao - bo)au 



2^ 



-ci + 



Cf. 



(31) 
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Finally, summing Eq.(pq) and Eq.(pll) one obtains 



4o7r 



2 ' 3a, 



(y. 

+ 2(oo - 6o)aie ( 1 + ^(3 + £2^(2))^^ 



47r 

(32) 



Eq.(32) is the final result. The 0{a^-) correction, represented by the term aie in Eq.(32) 



gets modified into {Cp = 4/3) 

aie ^ aie{l -^ '^ ^ ^ j . (33) 
Remembering that aie is written as (G^m^)/(\/2 87r^) in the OS and a/(47rs^)?n,^/(4m^) in 



the MS formulation, we see that one can take into account the leading QCD effect in and 



^2 



kh, or correspondingly and A;^, just by replacing in the one-loop 0{a-^) term mf by 



2 2 ( 1 asVr2-3\ 
"^f ^ I 1 - — 3 I • (34) 



Eq.(p3) coincides with the result found by Fleischer et al |12 



4 Conclusions 

In the previous Sections we have shown that the use of current correlation functions and their 
associated current algebra provides a very powerful and compact framework to discuss the 

2 _ 

0(a^7-) corrections in the vertex. 

2 

It is well known that the 0{a^j-) term should be finite by itself at the one-loop level. 
In fact there is no counterterm available to cancel any divergent contribution proportional to 
mf. The formalism of current correlation functions allows one to combine several Feynman 
diagrams and easily enforce the Ward identity that guarantee the finiteness of this term. The 

2 

actual calculation of the leading 0(a^-) correction becomes then trivial. 
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Furthermore the current algebra framework is very suitable to discuss the QCD corrections 

2 

to the one-loop 0(q— ^) term. In fact, because the strong and weak currents commute, the 
one-loop Ward identity in the electroweak sector is preserved. This is crucial for the structure 
of the divergent terms. Indeed no poles more divergent than 1/e are found (Cfr. Appendix B) 
because all the divergent terms are due to the QCD substructure. In comparison the standard 
two-loop calculation presents 0(l/e^), or even 0{\/e^) contributions if the IR divergencies 
are regularized using dimensional regularization [|l^ . From a practical point of view this fact 
greatly helps in performing the calculation. Especially the structure of the IR divergent terms 
is very simplified. We find only two different IR contributions (Fig. 4(a) and 4(b)) and their 
evaluation in the IR part is very similar to the computation of the one-loop QCD corrections 
to the Z66-vertex. 

2 

Concerning the physical significance of the 0(aas^-) calculation, we find that the QCD 
corrections to 0(a— ^) term in T{Z h + b) have opposite sign to that of the one-loop 
contribution. A similar situation happens in the electroweak parameters 5p and Ar. A 
detailed analysis shows that the inclusion of the perturbative QCD higher order effects 
in the vacuum polarization functions increases the prediction for rrit derived from current 
measurements. The values obtained for mt including the 0{aas) corrections are larger than 
those obtained using only the electroweak calculation by an amount between 5 and 10 GeV 
for 90 < mt < 200 GeV. Using a.^ = 0.118 @, Eq.(|^ gives as correction factor for mtAA%, 
that is consistent with the result of Ref. Q. In fact, a comparison between the correction 
found for 5p, i.e -^^^ ~ -2.860^^ and the value we obtained, ~ -2.290^^, 

shows that the two results are numerically comparable. Both corrections are quite large, 
however, it should be pointed out that they are scheme dependent. In our calculation, as 
well as in Ref. [P, the OS mass definition for the renormalized top mass has been used. If we 
employ instead a MS definition for the top mass, rht{p = rht), we have that the counterterm 
6mt has no finite part, and therefore in the evaluation of the diagrams in Fig. 5 only the 
terms proportional to 5 should be retained (Cfr. Appendix B). This changes the coefficient 
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multiplying ^ from — ^^"'^ to i.e. from ~ —2.290 to ~ +0.377, a much smaller 

number. Similarly the use of rht as renormalized mass in 5p changes the numerical coefficient 
in front of ^ from — ^^'g^^ to — g^^ ^ i-^- from ~ —2.860 to ~ +0.193, again a much smaller 
number. 

Finally we want to stress that our calculation of 0{aas-^^) terms and the one presented in 



Ref. [12| are completely independent. Because in Ref. |]12[ the IR divergencies were regularized 
using dimensional regularization it is not possible to make any intermediate check, not even 
at the level of T{Z ^ b + b). It is a welcome fact that our result coincides with the one 
obtained in Ref. IT! 
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Appendix A 



We wish to show that Eqs.(16) and ( |18a| ) contain correctly the field renormalization of the 
external quarks. The wave function renormalization constant Z;, for a massless b quark is 
equal to 

where Ilf, = ^bi^)/^\j^^Q with T,h equal to —i times the self-energy of the b. Expanding Zb 
up to 0{aas) considering only the contributions due to the $ and the gluon we have 

= 1 + + Hfe^ + • • • ~ 1 + ng' + nw + n(2) + 2ng'n« + • • • (A2) 
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where 11^' is the one- loop contribution due to the H^^^ the corresponding one for the gluon 
and n*^' the two-loop IPI mixed <I>-gluon term. Eq.(|A2|) tells us that any IPI self-energy 
diagram inserted in an external leg should be multiplied by a factor 1/2 while the total 
counting of the product of one-loop objects should be equal to 2. 

We consider first the vertex function V^g. Among the various contributions of the T*- 
product in Eq.( 15b| ) there are diagrams in which neither the strong currents nor the operators 



5 and enclose the vertex of the current. We begin by discussing these diagrams for 
the case of the b leg. They are depicted in Fig.6(a-g). The wiggly line in the figure stands 
for the q momentum absorbed by the operator. Again, as in Fig. 4 and Fig. 5, every graph 
represents two contributions with different momentum insertions. 

Indicating with Hg^^'o,) (a =a, b, c,d) the self-energy graph, divided by the relevant mo- 
mentum, inserted in the external leg of the diagram obtained by Fig. 6(a) attaching together 
the dashed lines, we can write the sum of the two diagrams represented by 6(a) as 



6(a) = lim -^f ^(^'Ot-^^^^^^ fej^ + )(/ - 4)] v{p) (A3) 



Remembering that H'^'has the form 11^^' = za^ and observing that the first term in the 
square bracket in Eq. ([A3D is zero we have 

6(a) = -if^^6(a)^(p07^ °°^/°^' a_z;(p). (A4) 

The factor l/2 2:g(a) is the correct contribution of the self-energy diagram associated with 
Fig. 6(a) to the field renormalization of the h. 

Let's consider now the graph 6(b). Each of the two diagrams of this kind, called 6(bi) 
and 6(b2), can be written as 

6(bi) = lim i^u(y)7^ "° + ^01^ ^-^1 [ng(^^j(p,p', ^)^ + ng(t,^)(p,/, g)(|i' - ^)] v{p) (A5) 



16 



where i = 1,2 and the l^e(b ) related to Hg^j^^ through 

^6(1) = ^ {Kb.M p'^ + Kb.) p'^ 9")) • (A6) 

It is clear that only the ng(^.j(p,p', g) term contributes in Eq. ([A5D . However, an explicit 
calculation shows that there is a relation between the Hg^^ j's, i.e. 

n6(bi) = ^6(1,2); Hgjj^^j = ng(,^^j. (A7) 

Therefore when we sum 6(bi) and 6(b2) and we use the relations (|A7|) and Eq.(|A^) we get 

6(b) = -i^^^z,^^^u{p')Y^^^p^a.v{p) (A8) 

that gives again the correct contribution to the field renormalization. 

The sum of the diagrams 6(c) and 6(d) can be discussed in a similar fashion as 6(b) 
obtaining 

6(c) + 6(d) = Qz6(c) + ^^6(d)) u{p')r^^^^a.v{p)- (A9) 

Therefore we conclude that V^^ takes correctly into account the two-loop IPI contribution 
to the wave function renormalization. 

Beside the IPI contribution Eg. ( |A2D contains the product of one-loop terms. Diagrams 
6(e-g) represent the corresponding contributions in V^^ for the b leg. Due to the limiting 
procedure we have that each diagram carries a factor 1/2. Therefore the sum (6e-g) plus 
the corresponding contribution of the b leg gives 3 times 11^' ^'g ^- There are other reducible 
contributions to V^g that do not give the correct renormalization factor. In particular dia- 
grams 6(h) is counted 1 in V^g because the limiting procedure does not affect it. We know, 
however, that it should be counted 1/2 because of the wave function renormalization of the 
external b leg. The same happens for the symmetric diagram where the gluon is on the 
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external b line. Therefore we have that in V^^ the product of one-loop $ diagrams times 11^^' 
is counted one time more than the correct contribution. The second term in Eq.(^) exactly 
cures this problem. The vertex function U^^ can be discussed similarly to diagram 6(h). The 
subtraction indicated in Eq.(|18a|) gives again the correct counting. 



Appendix B 

In this Appendix we list the contributions of Fig. 4 to and of Fig. 5 to c.t.. Defining 



2 m? 
5 = --27 + 21n47r-21n^ 
e 



and omitting the common factor 



4 (167r2)2 



we have 



^ 7r2 n2 \2 

4(a) = -^ + ^-ln^-ln2^-41n 



2 3 mf 

1 , A2 

o + In — 

2 mi 



4(b) = 4(c) = -5 - - + In 

''t 

4(d) = 4(5 + 13-6^(2) 

4(e) = 4(f) = -|<5-^ + 3^(2) 

4(g) = 4(h) = 45 + ^-^^(2) 

4(i) = 4(1) = +\- ^-m 

5(a) = 5(b) = 5(c) = -3(5 - 7. 

3 1 3 

5(d) = 5(e) = ^5 + - 
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The left-hand-sides in the above equations indicate the appropriate diagrams in Fig.4 and 
Fig. 5. In the same normalization the contribution of — n^^'g|-2?$ _ is equal to —4(b). 



q=q 
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Figure Caption 

Fig.l One-loop vertex diagrams contributing to the 0(a^^) corrections to the Z ^ b + b 
process. The figures schematically represent the sum of Feynman diagrams in which 
the propagators are attached in all possible ways to the external lines. 

Fig. 2 Momentum insertions associated with P$. The dotted line indicates a q momentum 
absorbed by the S2 operator. The dashed lines stand for the momenta absorbed or 
emitted by the and S operators. 

2 _ 

Fig. 3 Two-loop vertex diagrams contributing to the 0(aas^^) corrections to the Z ^ b + b 
process. The meaning of the circles is explained in the text at the beginning of Section 
3. 

Fig.4 Types of graphs contributing to V^g. Every graph represents two diagrams with 
different momentum insertions (see text). The meaning of the lines is as in Fig.2. 

Fig.5 Types of diagrams belonging to the counterterm contributions {c.t.). Every graph has 
the same meaning as in Fig.4. 

Fig.6 Subset of graphs contributing to V^g- The meaning of the graphs is the same as 
in Fig.4 and Fig.5. The wiggly line stands for the q momentum absorbed by the 
operator. 
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